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Abstract 

A so called ^-umbral extensions of the Stirling numbers of the second 
kind are considered and the resulting Dobinski-like various formulas 
- including new ones - are presented. These extensions naturally en- 
compass the two well known q-extensions. The further consecutive 
ip- umbral extensions of Carlitz-Gould-Milnc g-Stirling numbers are 
therefore realized here in a two-fold way. The fact that the umbral 
(/-extended Dobinski formula may also be interpreted as the average of 
powers of random variable X q with the g-Poisson distribution singles 
out the g-extensions which appear to be a kind of "singular point" in 
the domain of -0-umbral extensions as expressed by Observations 2.1 
and 2.2. Other relevant possibilities are tackled with the paper's closing 
down questions and suggestions with respect to other already existing 
extensions while a brief limited survey of these other type extensions is 
being delivered. There the Newton interpolation formula and divided 
differences appear helpful and inevitable along with umbra symbolic 
language in describing properties of general exponential polynomials of 
Touchard and their possible generalizations. Exponential structures or 
algebraically equivalent prefabs with their exponential formula appear 
to be also naturally relevant. 

MCS numbers: 05A40, 11B73, 81S99 

Key words: extended umbral calculus, Dobinski type formulas, Graves- 
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1. In the g-extensions realm 

At first let us make a remark on notation (see also Appendix) . t/j is a number 
or functions' sequence - sequence of functions of a parameter q. tp denotes an 
extension of (^r)n>o sequence to quite arbitrary one (the so called - "admissi- 
ble" [1, 2]). The specific choices are for example : Fibonomialy-extended se- 
quence (-p3-|)n>o {{Fn) - Fibonacci sequence ) or just "the usual" V-sequence 
(^r) n >o or Gauss g-extended (^r)n>o admissible sequence of extended um- 

bral operator calculus, where n q = \^ q and n q \ = n q (n — l) g !, q \ = 1 - see 
more below. With such type extension we may " ip- mnemonic" repeat with 
exactly the same simplicity this what was done by Rota forty one years ago. 
Namely forty one years ago Gian-Carlo Rota [3| proved that the exponential 
generating function for Bell numbers B n is of the form 

oo n 

Y, ^jB n = exp(e* - 1) (1) 
n! 

ra=0 

using the linear functional L such that 

L(X*) = 1, n > 0. (2) 

Bell numbers (see: formula (4)in [3]) are then defined by 

L(X n ) =B n , n> 0. (3) 

The above formula is exactly the Dobinski formula [1] if L is interpreted 
as the average value functional for the random variable X with the Poisson 
distribution where L(X) = 1. As a matter of fact it is Blissard calculus 
inspired umbral formula [H] (see [5] for umbral nature of Poisson random 
variables and the introduction in [5] for historical remarks on Blissard's 
calculus roots). On this occasion let us recall that the Stirling numbers of 
the second kind are relatives of the Poisson distribution in the known way. 
Namely if A is a random variable with a Poisson distribution with expected 
value A, then its n — th moment is the exponential polynomial <p n value at 
A i.e. 

n 

E{X n ) = ^ n {X)=Y J {lY k - 

k=0 

Hence in particular, the n — th moment of the Poisson distribution with 
expected value 1 is precisely the number of partitions of a set of size n i.e. 
it is the n — th Bell number (this fact is Dobinski's formula as stated by 
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the formula (3)). The formula (3) is tempting to be ^-extended somehow 
as the ?/>-Poisson process distribution is known [2, 1]. Before doing this let 
us remind that recently an interest to extensions of Stirling numbers and 
consequently to Bell numbers was revived among "(/-coherent states physi- 
cists" H3 |H] with several important generalizations already at hand such 
as in [9, 10, 11]. The merit of such applications is in that the expectation 
value with respect to coherent state (7 > with | -y | = 1 of the n-th power 
of the number of quanta operator |7| |H] is "just" the n-th Bell number 
B n and the explicit formula for this expectation number of quanta is "just" 
Dobinski formula B|. The same holds for (/-coherent states case |B| i.e. the 
expectation value with respect to (/-coherent state [7 > with [7] = 1 of the 
n-th power of the number operator is the n-th (/-Bell number [HJ |H| defined 

as the sum of (/-Stirling numbers j ^ j introduced by Carlitz and Gould and 

recently exploited among others in j!21 |SJ [3 El- Note there then that for 
the two standard (/-extensions of the Stirling numbers of the second 
kind we have as the first ones the (/-Stirling numbers: 



k=0 H 

where x q = ^jzh an d %q = x q (x — l) q ...(x — k + l) q and then the second ones 
called (/"""-Stirling numbers. Both correspond to the tp sequence choice in the 
(/-Gauss form (— j) n >o- Here the (/^"-Stirling numbers of the second kind are 
introduced as coefficients in the famous Newton interpolation formula (Liber 
III, Lemma V, pp. 481-483 in [13]) now applied to the polynomial sequence 
( e «-)n>o> e "( x ) = x n ,n> 0, i.e. 

n 

x n = ^Z{ n k } Xk(x), i.e. {^} = [0,l q ,2 q ,...,k q ;e n ], (5) 

k=0 ' ' 

where Xk(x) = x(x — l q )(x — 2 q )...(x — [k — l] q ), and 

[x ,x 1 ,X2,...,x k -i;f]-[xi,X2,...,x k ;f] 

[x Q ,x 1 ,X2,...,x k ;j\ = 

xq - x k 

denotes the k — th divided difference with 

[xo,x\U\ = • 

XQ - Xl 

These two kinds of now classical g-extensions of Stirling numbers of the 
second kind as defined by (4) and (5) are related in a simple way through 
re-scaling [14]. They satisfy the known respective recurrences: 
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H 1=0 V 7 9 y 

From the above it follows immediately that corresponding g-extensions of 
B n Bell numbers satisfy respective recurrences on their own: 

B q {n + l) = j^( n \ q l B g (l);n>0, 
i=o V A 

where 

fc=0 y 

Different definitions via (4) and (5) equations correspond consequently to 
different (/-counting [Tl]. With any other choice out of countless choices of 
the i/j sequence the equation (5) becomes just the definition of ^-Stirling 
(vide " Fibonomial-Stirling" ) numbers of the second kind | ^ | and then 

■i/i^-Bell numbers B~(ip) are defined as usual as sums of the correspond- 
ing Stirling- like numbers - where now Xk( x ) m (5) is to be replaced by 
?Pk(x) = x(x — l^)(x — 2, l p)...(x — [k — 1]^)- These -(/^-Stirling numbers of the 
second kind for q case identified as Comtet numbers in Wagner's terminol- 
ogy H4j satisfy familiar recursion and are given by familiar formulas to 
be presented soon. The extension of definition (4) of the q- Stirling numbers 
of the second kind beyond this (/-case i.e. beyond the ip = (^ri}n>o choice 
is not that mnemonic at all and the problem of recursion appears. This 
part of alternative treatment is to be considered later on after we exploit a 
little bit more some consequences of (4). Namely - due to (4) one imme- 
diately notices that the expectation value with respect to (/-coherent state 
|7 > with I7I = 1 of the n-th power of the number operator is exactly the 
popular (/-Dobinski formula which can be given Blissard calculus inspired 
umbral form - like in (3). It is enough to apply to (4) L q - the average value 
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functional for g-Poisson distribution [1, 2]. The formula thus obtained may 
be also treated as a definition of g-extended Bell numbers B n (q) 

L q {X n q ) =B n (q), n>0 (6) 

due to the fact that this linear functional L q interpreted as the average value 
functional for the random variable X q with the g-Poisson distribution [1] ( 
L q (X q ) = 1 ) satisfies 

L q {Xf) = 1, n > 0. (7) 
Then with the g-exponential polynomials 

n 

(fn(x,q) = ^2 {^} s * (q - exp - pol - I) 

one obtains for x = 1 in correspondence with L q (X q ) = 1 the q - formula of 
Dobinski type (compare with [12] and see (5.28) in [15]): 

ip n (l,q)=B n (q)=L q {X n q )=e q l Y.T-V n ^ °' e ? = ^Wl" 1 - 

k=o K o- 

We arrive to this simple conclusion using Jackson derivative difference op- 
erator in place of D = d/dx in q =1 case and the power series generating 
function G(t) for g-Poisson probability distribution: 

\n 

p n = [expgX]- 1 — ,G(t) = £>„i n , (8) 

H 1- „>0 

where exp q X = exp^X = Y, n >o ^n(q)t n , ipn(q) = ^t- Naturally 
d n G(t) 

Pn = [^p]t=o, [d q G(t)] t=1 = 1 for X = l. (9) 

In order to arrive at the g-Dobinski formula apply (7) to (4) with (6) in mind. 
(As for f/'-Poisson probability distribution - see [1,2].) There are many ways 
leading to g-extended Stirling numbers according to their weighted counting 
interpretation i.e. various statistics are counted by g-Stirling numbers of the 
second kind F° r example w(ir) = q cross ( n )^ w(ir) = from |17j 

give the Carlitz-Gould- Milne g-Stirling numbers | & | " after being summed 
over the set of A:-block partitions while w(ir) = q mn ( w ) from 18^ gives rise 
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to the Carlitz- Gould -Milne (/^-Stirling numbers | £ j after being summed 

over the set of fc-block partitions (see also maj and maj'~ statistics in ^5] 
as well as other statistics [14] giving also both extensions). The choice of 
weig ht w(n) = g*W HO] gives rise to another statistics counted by new kind 
of (/-extended Stirling numbers of the second kind. Let us consider - for 
the sake of illustration this example from [201 m more detail. This is the 
example of weighted counting ^o-en n °f partitions of [n]. The weight 
w of such a partition ir 6 II n is determined by inversions' i function in the 
form w(tt) = q l ^\ Here Tl n denotes the lattice of all partitions of the set [n] 
while A Ut k denotes the family of all k -block partitions. A fc-block partition 
7T G A n ^ C n n is represented in the standard form: ir = B\j 'B2/ '.../ 'B^ with 
the convention that maxB\ < maxBi < ...,maxBk = n. For i € [n] let 
hi denotes a number of a block to which i pertains. Define an inversion 
of partition ir to be a pair < i,j > such that bi < bj and i > j. The 
inversion set of tt is I(ir) = {< i,j >; < i,j > is an inversion of ir}. Then 
i(ir) = [/(vr)| and the inversion (/-Bell numbers are naturally defined as 

bt (?) = E ^ = E E a 400 

while inversion (/-Stirling numbers of the second kind are identified with 

The inversion (/-Bell number B™ v (q) is the generating function for I(s) = 
the number of all partitions from Tl n with s of the above inversions because 

B l n nv (q) = E 9* W =E^ E 1 = E W- 

crGilrj s>0 TTeU n ,i(n)=s s>0 

Recursions for both inversion (/-Bell numbers and inversion (/-Stirling num- 
bers of the second kind are not difficult to be derived. Also in a natural 
way the inversion (/-Stirling numbers of the second kind from [20] satisfy 
a (/-analogue of the standard recursion for Stirling numbers of the second 
kind to be written via mnemonic adding "g" subscript to the binomial and 
second kind Stirling symbols in the the standard recursion formula i.e. 

rn + li™ ^ fn\ (n-l^i™ 
1=0 v 7 1 
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Another g-extended Stirling numbers much different from Carlitz "g-ones" 
were introduced in the reference |5J from where one infers j^2] the cigl- 
analog of (6). Let IT denotes the lattice of all partitions of the set {0, 1, .., n — 
1}. Let 7r £ LT be represented by blocks ir = {B Q , B±, ...Bi, ...}, where B Q is 
the block containing zero: 6 B a . The weight adapted by Cigler defines 
weighted partitions' counting according to the content of B Q . Namely 

Therefore ^ttgii q c%al ^ = B n {q). Here ^4 nj fc stays for subfamily of all fc-block 
partitions. With the above relations one has defined the cigZ-g-Stirling and 
the cigl-q-TSe\l numbers. The cig/-g-Stirling numbers of the second kind are 
expressed in terms of g-binomial coefficients and q = 1 Stirling numbers of 
the second kind |2J as follows 

r n + 1 1 ci 9 l v-^ / n\ (n-t+i\ r n — 1 1 ci 9' 

{ t }, =g(J,' ( 1 VJ. i-so.*si- 

As seen above these are new q-extended Stirling numbers. The correspond- 
ing cig/-(/-Bell numbers recently have been equivalently defined via cigl-q- 
Dobinski formula |22| - which now in more adequate notation reads: 



L{X qn ) = B n (q), n>0,X q " =X(X + q-l)...(X -l + q 11 - 1 ). 

The above ci(//-g-Dobinski formula is interpreted as the average of this spe- 
cific n — th cigl-q-power random variable X q " with the q = 1 Poisson dis- 
tribution such that L(X) = 1. For that to see use the identity by Cigler 



x(x-l + q)...(x - 1 + g"" 1 ) = { n k Y^ xk 

k=0 



2. Beyond the g-extensions realm 

The further consecutive ^-umbral extension of Carlitz-Gould g-Stirling num- 
bers j ^ j and |^| is realized two-fold way - one of which leads to a 
surprise in contrary to the other way. 
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2.1. The first way 



The first "easy way" consists in almost mnemonic sometimes replacement 
of q subscript by ifi after having realized that in equation (5) we are dealing 
with the specific case of the so called Comtet numbers |141ll5j (Comtet L. in 
N ombres de Stirling generaux et fonctions symtriques C.R. Acad. Sci. Paris, 
Series A, 275 (1972): 747-750 formula (2) refers to Wronski). This array of 

Stirling-like numbers |^| " "ctlephs de Wronski" as Comtet refers to it or 

these Comtet numbers in terminology of Wagner [14, 15] or as a matter of 
fact [13] these Newton interpolation coefficients for e n ,n > i.e. divided 
differences [0, 1^, 2^, k^; e n ] are defined accordingly as such coefficients - 
below. 

n 

x n = Y,{l} ~M x )> n ^°> ( 10 ) 

fc=0 

i.e. equivalently (recall that e n (x) = x n ,n > 0) 



(?) = [0, V, 2^, k^;e n ] = Y) ? n l fj , n > 0, (Newton) 
where 

ipk(x) = x(x - h/,)(x - 2^)...(x -[k- %) 

and ij) s denotes the first derivative. Let then / = (f n )n>o be an arbitrary 
sequence of polynomials. In the following we shall call S(f; n, k) defined 
below 

[do,d 1 ,d 2 ,...,d k ;f n ] = S(f){di}i> ,n,k) (N-W-C Stirling) 

the Newton- Wronski-Comptet Stirling numbers (N-W-C for short)- compare 
with Appendix A. 2. 

The - Stirling numbers "f & } defined by (10) are specification of N—W— 

C Stirling array for which we naturally define ^-exponential polynomials 
(p n (x,ifi) as follows 

n 

<p~(x,ip) = ^[0, l^,2^,...,k^;e n ]x k , n > 0. (ip~ - exp - pol) 

k=0 

Note the trivial but important fact that in the N-W-C Stirling numbers 
case we are dealing with not equidistant nodes' interpolation in general and 
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note that (Rescal) from the subsection 2.2. below is no more valid beyond 
(/-extension case - both with an impact on the way to find out the Dobinski- 
like formulae - see more below. 

As a consequence of (10) we have "for granted" the following extensions of 
recurrences for Stirling numbers of the second kind: 

where |g| = 6 n> o, |^| =0, k > n; and the recurrence for 
ordinary generating function reads 



G k^) = T ^ V G^_ 1 (x), k>l, (12) 



where naturally 



n>0 

from where one infers that 



* (1 - l^x)(l - 2^x)..\l - k^x) 

Hence we arrive in the standard extended text-book way [221 a ^ t ne following 
explicit new formula (compare with (2.3) in [15]) 



v r= i \ 



n>k>0, 



(14) 

where 



£(-l)*-(J*)r$; n,^>0 



r=l 

is readily recognized as the ^-extension of the formula for surjections in its 
- after inclusion-exclusion principle had been applied - form. 

Expanding the right hand side of (13) results in another explicit formula 
for these ■i/'-case Newton- Wronski-Comtet array of Stirling numbers of the 
second kind i.e. we have 

{fc} = (*i)^(*2)v>-(*n-fc)i/>; n>k>0 (15) 

l<ii<i2<...<i„_fc<£! 
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or equivalently (compare with [13, 14]) 

n i ~ 
k Jip 



" V - E lf2f...^; n>k>0. (16) 



di+d2+...+df c =n— k, di>0 

N-W-C case ip~ - Stirling numbers of the second kind being defined equiv- 
alently by (10), (Newton), (14), (15) or (16) yield N-W-C case ^~ - Bell 
numbers 



n n 

fc=0 fc=o 

Naturally 3! functional L~ such that on the basis of persistent root poly- 
nomials ipk(x) it takes the value 1: 

L~(Vfc(s)) = 1, & > 0. 
Then from (10) we get an analog of (3) 

Problem: which distribution the functional LT is related to is an open tech- 
nical question by now. More - the recurrence for B£(ip) is already quite 
involved and complicated for the g-extension case (see: the first section)- 
and no acceptable readable form of recurrence for the ^-extension case is 
known to us by now. 

Nevertheless after adapting the standard text-book method |23| we have 
the following formulae for two variable ordinary generating function for 
| & | Stirling numbers of the second kind and the ^-exponential gener- 
ating function for B~(ifi) Bell numbers 

c£(*,y) = X>~0M* n , (17) 

n>0 

where the ip- exponential polynomials (p~(i(),y) 

n 

=£{!}>' 

fc=0 

do satisfy the recurrence (compare with formulas (28) in Touchard's [24] 
from 1956) 
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VnW.I/) = [!/( 1 + ^]CiW.!/) n>\, 

hence 

<p~(i>,y) = [y(l + d^} n l, n>0. 

The linear operator cL, acting on the algebra of formal power series is being 
called (see: [1, 2] and references therein) the " ^-derivative" as d^y n = 

The - exponential generating function 

B^x) = Y j B:^)- 1 W~-e.g.f.) 

n>0 

for -B~(V0 Bell numbers - after cautious adaptation of the method from 
the Wilf's generatingfunctionology book [23] can be seen to be given by the 
following new formula 

W=£e(^,r)^M (18) 

r>0 

where (see: [1,2] and references therein) 

n>0 Y 

while 

^• r) =£^b <19) 

and the new Dobinski - like formula for the ^-extensions here now reads 

The ^"""-exponential polynomials are therefore given correspondingly by 
9?~(V', x) = e (V') r)—^x r . — exp — pol — II) 



r>0 



In the case of Gauss (/-extended choice of (^-f)n>o admissible sequence of 
extended umbral operator calculus equations (19) and (20) take the form 



k=r 
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and the new N-W-C case q^-Dobinski formula is given by 

= X>(ff.Orf, ( 22 ) 



r>0 q ' 



which for q = 1 becomes the Dobinski formula from 1887 [4]. Note the 
appearance of re-scaling factor q~^ in (21). In its absence we would get 
not q^-Dobinski but q-Dobinski formula 

B n (q) = ^j r J2^ i (q- Dobinski) 

- see [15] and formula (5.28) there coinciding with N-W-C case of Dobinski 
formula after re-scaling in correspondence with (Rescal) below in subsection 
2.2. Correspondingly we would the have not (qr~ — exp — pol) formula but 
(q — exp — pol) formula: 



n n 

<Pn(x,q) = ^<?( 2 )[0, l q ,2 q , k q ;e n ]x k = ^ {"} x k . (q - exp - pol) 

k=0 k=0 q 

The interpretation problem. Combinatorial interpretations of the known 
up to now various ^-extensions of Stirling numbers of both kinds - are briefly 
reported on in the Appendix. The problem of how eventually one might 
interpret - beyond the q-extensions' realm - for example the ^"""-Dobinski 
formulae (20) and (22) also in the Rota-like way represented here by equa- 
tion (3) we leave opened - see the discussion in Appendix A. 2 and A. 3. II. 
Naturally there exist a unique linear functional such that 

L^ k (x) = x(x - l q )(x - 2 q )...(x -[k- %)) = 1, k > 0. 

It is also to be noted that in the exceptional case of ^-extensions and only 
for ^-extension we have equivalence of | \ | and | \ | by re-scaling. 

For the latter ^-Stirling numbers we have Dobinski formula and simulta- 
neously g-Poisson average functional interpretation as represented by the 
definition (6). Namely - recall the fact that the linear functional L q there is 
interpreted as the average value functional for the random variable X q with 
the g-Poisson distribution which is specific case of the ?/>-Poisson distribu- 
tions from [2, 1]. Note again that this re-scaling takes place for (/-extensions 
and only for g-extension. This is so because the relation 
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"V - Aty = f(k)(n - k)^, 1^ = 1 

holds for and only for f(l) = q when it becomes the identity (2.6) from 
[25] i.e. f(k) = q k . It is our conviction that this is the very reason that q- 
extensions seem to appear as a kind of "a bifurcation point" in the domain of 
V>-umbral extensions. This conviction is supported by the corresponding 
considerations in [26]- section 3 - on possibility of -^-analogue of the so 
called "quantum" g-plane formulation of g-umbral calculus. The parallel 



treatment of the Newton- Wronski-Comtet 



n 
k 



Stirling numbers of the 



first kind is now not difficult (consult [25, 12, 8] for example and Wagner's 
recent treatment of the well established (/-case in [14]). 

In the inversion-dual way to our equation (10) above we define the tp~- 
Stirling numbers of the first kind as coefficients in the following expansion 



r=0 



k 



(23) 



where - recall ipk(x) = x(x — l 1 p)(x — 2^)...(x — [k — 1]^,). (Attention: see 
equations (10)- (16) in [8] and note the difference with the present definition). 
Therefore from the above we infer that 



E 



=o 



k 



(24) 



Another natural counterpart to ^"-Stirling numbers of the second are 
Stirling numbers of the first kind defined here down as coefficients in the 
following expansion (" c" because of cycles in non-extended case) 



r=0 



k 
r 



(25) 



where - now = x ( x + ^4>)i x + 2^,). ..(a; + [k — 1]^)- These are to be 

studied elsewhere. 

On interpretation. For possible unified combinatorial interpretations of 
binomial coefficients of both kinds, the Stirling numbers of both kinds and 
the Gaussian coefficients of the first and second kind - i.e for the specific 
choices of ip = (^)n>o - see [27, 28]. As for q-analogue of Stirling cycle 
numbers see [29] and Sect. 5.3. in [30]. The problem of eventual combina- 
torial interpretation of other -^-extensions (vide Fibonomial - for example) 
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- remains opened. 



2.2. The second way. 

We shall come over now to inspect the outcomes of the second way - 
announced at the start of this section - after having realized that in the 
equation (4) we are not dealing straightforwardly with Newton- Wronski- 
Comtet array of Stirling-like numbers ^3 El El - except for the q = 1 case 
"extension"- of course. Though it is to be noted that still the re-scaling 
takes place for g-extensions and only for (/-extension. Not beyond. (This 
is so because the relation 

n<p - kjf, = f(k)(n - k)^, 1^ = 1 

holds for and only for /(l) = q when it becomes the identity (2.6) from 
Gould's [25] i.e. f(k) = q k ). Thus after the above Gould re-scaling we 
would recover Newton-Wronski-Comtet array of Stirling-like numbers re- 
scaled - anyhow i.e. 

n 

y n = E«~®{fc} y "' {Rescal) 

k=o q 

where y = x q = and ^fc(y) = y(y - l q )(y - 2 q )...(y - [k - l] q ). At first 
let us recall that the definition (4) of g-Stirling numbers of the second kind 
| ^ \ is equivalent to the definition by the recursion 

where {™}^ = <5 ni0 , {fc} g = °> k > n - 

These in turn is equivalent to (just use the standard Q-Leibniz rule [1, 2, 
31] for Jackson derivative d q ) 

(xd q r=j2{ n k } (27) 

fe=o q 



where 




= 0, k > n. Here x denotes the multiplica- 



tion by the argument of a function. The formula (27) is a special case of the 
typical for GHW algebra [32, 1, 2, 33] expression investigated by Carlitz in 
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1932 [34] (compare with GHW formulae (1), (31), (32) in [7] and see also 
[35]). The idea now is to extend eventually the definition by equation (27) 

via replacing (/-extended operators by the corresponding -(/'-extended ele- 
ments of the Graves-Heisenberg-Weyl (GHW) algebra representation [32, 
1, 2, 33]. However we note at once (see: Appendix for particulars of the 
up-side down notation) that the two following observations hold. 
Observation 2.1 The equivalent definitions (28) and (29) 

n 

(xd f r = j2{ n k }/ k9 i ( 28 ) 



fc=0 



where {o} ?/ , = ^ n > ' {fc}^, = ' k > n and 



±{114 (») 

k=0 



lead to the one first order recurrence of the type (26) for and only for 
(/-extension. 

This again is so because the relation 

rty - Aty = f(k)(n - k)^, 1.0 = 1 

holds for and only for f(l) i.e. for (/-extension, where it becomes the identity 
(2.6) from [25] i.e. f(k) = q k . The next observation now comes as a would 
be surprise. 

Observation 2.2 The equivalent definitions (28) and (29) have no non- 
trivial realizations beyond the ^-extension case. 

In order to arrive at this observation let us act appropriately on x N 
N > monomials by both sides of the GHW algebra representation defi- 
nition (28) thus getting an infinite sequence of recurrences 



k>0 v k>0 



k U v ^YkH 



with no nontrivial solutions as spectacularly evident with the choice of - 
for example - Fibonomialy-extended sequence (^j)n>o ((Fn) - Fibonacci 
sequence ) unless ip = (^-\)n>0- And for this and only for this choice 
V' = (^)n>o we have 

N q = q k ((N-k) q + k q 



15 



which after being applied in (30) results in one recurrence which is exactly 
the recurrence (26). 

As expected - the equation (29) becomes equivalent to the one first order 
recurrence of the type (26) for and only for (/-extension. 

Closing remark. We see that the Carlitz-Gould (/-Stirling numbers { & | 

make (/-umbral extension to appear as a kind of "a bifurcation point" in the 
domain of respective -0-umbral extensions. This in statu nascendi convic- 
tion is also supported by the corresponding considerations in [26]- section 3 
- considerations about possibility of 'i/'-analogue of the so called "quantum" 
(/-plane formulation of g-umbral calculus. As for eventual second way's ip- 
extensions beyond the (/-extension case where the rescaling does not take 
place we are left with an opened problem how to eventually find the way to 
get round this inspiring obstacle. The selective comparison of the presented 
umbral extensions of Stirling numbers, Bell numbers and Dobinski-like for- 
mulas with other existing extensions (as well as relevant information in brief) 
serves the purpose of seeking analogies and is to be find in the Appendix 
that follows now. 



Appendix - for remarks, discussion and brief comparative re- 
view of ideas. 

A.l. Notation. 

The necessary commutation relations' representation for the GHW (Graves- 
Heisenberg-Weyl) algebra generators is provided in [31, 1, 2, 33]. Applica- 
tions of these might be worthy of the further study [25, 36, 37]. The sim- 
plicity of the first steps to be done while identifying general properties of 
such ^-extensions consists in notation i.e. here - in writing objects of these 
extensions in mnemonic convenient upside down notation [J, [2] 



VV-l) i / -i \ i n ip(x-l) , . 

— = ru,, Thhl = 7hb{n- lUl, n > 0, x$ = — ~r\ — , (31) 

x \ = x i>( x - l )ip{ x - 1)rji—{x + (32) 

*,(* - 1V..(* - k + 1), = ^ (x _ lMKx _ fc + 1) - (33) 

If one writes the above in the form = ^fey^ = $(x) = <& x = x$, one sees 
that the name upside down notation is legitimate. You may consult [1, 2, 26, 
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33, 36, 37] for further development and usefulness of this notation. In this 
notation the ^-extension of binomial incidence coefficients read familiar: 

/n\ _ ?V _ / n \ 
\k) ip k^\{n — k)^,\ V n — k) ^ 

A. 2. Discussion, remarks, questions. 

g-umbral extensions are expected to be of distinguished character - also due 
to what was stated in Section 2. Being so they pay to us with simplicity 
of formulae and elegance of q - weighting combinatorial interpretations and 
thus various statistics of the combinatorial origin. 

Because of that and because of the importance of g-umbral extensions in 
coherent states mathematics we here adjoin a remark on simplicity based 
on the remark of Professor Cigler (in private) . 

Namely Katriel indicates in the very important source paper [Sj that his 
derivation of the Dobinski formula is the simplest. And really it is simple 
and wise. Possibly then this may be occasionally and profitably confronted 
with the also extremely simple derivation by Cigler (see p. 104 in [38]) 
based on GHW-algebra properties. Let then x denotes the multiplication 

by x operator while D denotes differentiation - both acting on the prehilbert 
space P of polynomials. Then due to the recursion for Stirling numbers of 
the second kind and the identity (operators act on P) 

i(D + 1) = ——(xD) exp (x) 

exp (x) 

one defines in GHW - algebra manner the exponential polynomials 

n 

n>0, ip n (x) = { n {ExPol) 

introduced by Acturialist J.F. Steffensen [39, 40] (see: Bell's "Exponen- 
tial polynomials" in umbra-symbolic language [41] p. 265 and his symbolic 
formula (4.7) for now Bell numbers). These exponential polynomials were 
substantially investigated by Touchard in Blissard umbra-symbolic language 
[24] . Here now comes the GHW-definition [38] of these basic polynomials 

<Pn(x) = l -r^{xD) n exp (x) 

exp (x) 
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resulting in the formula which becomes Dobinski one for x = 1 i.e. 



Note: The (/-case as well as ip-case formal mnemonic counterpart formulae 
are automatically arrived at with the mnemonic attaching of q or ip indices 
to nonnegative numbers |2j - vide: 



which for ip = (— } n >o an d x = 1 becomes the well known g-Dobinski 
formula as of course (p n (x = L?) = B n (q) - see in [15] the formula (5.28) 
and note that this is not g~-Dobinski formula (22) as noticed right after 
(22). As for eventual second way's ^-extensions beyond the g-extension case 
where the rescaling does not take place - we are left with an opened problem 
how to eventually find the way to get round this inspiring obstacle. Perhaps 
instead of the second beyond the g-extension way we might follow Alexander 
the great in his Gordian Knot problem solution and define S(ip, n, k) as 
follows (whenever one may prove that the object being defined is really a 
polynomial): 



An alternative good idea perhaps would be an attempt to V'-extend the cele- 
brated Newton interpolation formula ( use instead D, then exp^ instead 
of exp and then you will be faced with ^-Leibniz rule application problem 
though... see [1, 2, 33] for Leibnitz rules). Let us then make - also for the 
sake of comparison with existing knowledge - let us then make us wonder 
on the intrinsic presence and assistance of Newton interpolation which cor- 
responds to the first "easy" way as described in Subsection 2.1. 
The intrinsic presence and assistance of Newton interpolation for- 
mula in derivation of Dobinski formula for exponential polynomials and their 
binomial analogues was underlined and used in [42] for specific presentation 
of the q = 1 case from the umbral point of view of the classical finite opera- 
tor calculus. In [42] a Dobinski- like formula was derived being as a matter of 
fact the particular ("binomial") case of formula (30) from Touchard's 1956 





(34) 




(S(ip) — exp — pol) 



18 



year paper [24]. In more detail. Choosing any binomial polynomial sequence 
(b n ) n>0 consider its Newton interpolation formula 

n 

K(x) = 2}0,l,2,...,fc;& n ]a£. 

Then apply an umbral operator sending the binomial basis {x—) n>0 of delta 
operator A to the binomial basis (x n ) n >o of delta operator D. Then use 

[0, 1, 2, k; b n ] = = p ^k-l^ {NeWtOU ~ SUrlin9) 

so as to arrive (thanks to binomial convolution) at Dobinski like formula 
from [42] i.e. 

where ip is the umbral symbol satisfying [24] 

ip n+ i = x((p + l) n , = x k . (Touchard) 

In order to see that this is just the particular ("binomial") case of umbra- 
symbolic formula (30) from Touchard's 1956 year paper [24] just choose in 
Touchard formula (30) the arbitrary polynomial / to be any binomial one 
b n = f. Then f(<p) = b n (tp) = b n (ip(x)) is binomial also and we have 

f ^ = elplV) ^ bn ^. X " (Dobinski - Touchard) 

Equidistant nodes Newton's interpolation array of coefficients [0, 1, 2, k; b n ] 
- here the connection constants of the general exponential polynomial p n (x) = 
b n ((p(x)) are to be called in the following the Newton-Stirling numbers of 
the second kind and are consequently given by 

Pn(x) = — l -- £ ~rr~~ = Et°' L 2 ' -» fc ' (N-S-Dob) 
exp(z) ^ k\ ^ 

where (b n ) n >o is any sequence of polynomials. These are - in their turn - 
the special case of N-W-C Stirling numbers. 
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Coherent States' Example I. Take the b m (x) = f(x) in the (Dobinski- 
Touchard) formula to be of the form resulting from normal ordering problem 
(see A.3.II. - below) i.e. let (see: [10]) 

n 

f{x) = b ns (x;r,s) = ]J[x + (j - l)(r - s)] a 
i=i 

Then we get (2.8) from [10] i.e. 

-» k 7 

[0, 1, 2, k- b ns (...; r, s)] = - ]T (-1)*-'M*; r, s) ( f ) = 5 r , s (n, fc) 

i=s 

becomes the definition of the generalized Stirling numbers (see A.3.II. - 
below), which appear to be special case of general Newton-Stirling num- 
bers of the second kind. (Here b ns (...;r, s)(x) = b ns (x;r, s).) Naturally 
the Dobinski-like formula (2.1) from [10] for exponential polynomials de- 
termined by [0, 1, 2, k; b ns (.; r, s] = S r ^ s (n,k) is special case of (N-S-Dob) 
Dobinski-like formula with counting adapted to the choice / = b ns . Along 
with Bell numbers' sequence or Bessel numbers's sequence this special case 
of Newton-Bell numbers' sequence 

ns 

B r ,s(n) = y^S rtS (n,k) 

l=s 

is a moment sequence [43]. 

Example II The next example of Newton-Stirling numbers d n k comes from 
the paper [44] on interpolation series related to the Abel-Goncharov problem. 
There the divided difference functional is applied to e n yielding d n ^ 
accordingly: 

12 k - 1 

£±k&n — IVj ~r, jr , lj Gn\ — 0-n,k- 

The general rules for Newton-Stirling arrays allow us to notice that 

< k = [0 ' v v nr ' 1; en] = ¥ 5> 1)fc ~ r £ ; n - k - °' 

r=0 V y 

hence for corresponding exponential polynomials we have 

n , j. k 



fc=0 r=0 



20 



in accordance with the fact [44] that k n k d n ^ = |& j- Derivation of the 
Dobinski-like formula we leave as an exercise. 

On ^-extension. A ■(/'-extension of the above Touchard's symbolic defini- 
tion of exponential polynomials would start with the defining formula 

<p n+1 = x{ip+^ l) n , tf£ = x k . - Exp - Pol) 

resulting in analogous umbra-symbolic identities and with corresponding 
Dobinski-like formula as (35) below, where b n = e n . Compare these with 
(10) from where we have for this case of b n {x) = e n {x) = x n ,n > the 
Newton interpolation formula 

n 

x n = ^[0, 1^, 2ip, k^;e n ]ipk{x), n > 0. 

k=0 

For the meaning of the ^-shift "+/ see [1, 2, 26, 31, 33]. This we shall 
develop elsewhere. Meanwhile let us continue the limited review of other 
extensions. 

Plethystic Stirling numbers' extension The above umbral extensions 
as well as the other extensions to be mentioned in what follows are to be 
confronted with inventions of plethystic exponential polynomials, plethystic 
Stirling numbers of the second kind and plethystic extension of Bell numbers 
from [45] which constitutes an advanced and profound way to reach the merit 
of the finite operator calculus representations - this time realized with vector 
space of polynomials in the infinite sequence of variables. In [45] Mendez 
had derived profits from Nava's combinatorics of plethysm then developed 
by Chen to become an elegant plethystic representation of umbral calculus 
so as to find out also umbral inverses of plethystic exponential polynomials 
and related plethystic Stirling numbers of the first kind (for references see: 
[45]). The plethystic exponential polynomials are then there expressed via 
Dobinski-like (plethystic ) formula ( see: (35) in [45]) and the plethystic 
Stirling numbers of the second kind (see: (38) in [45]) are expressed via for- 
mula extending the formula for Stirling numbers of the second kind resulting 
from the formula for surjections in its after inclusion-exclusion principle had 
been applied form. Whether -0-extension of plethystic constructs as above is 
interesting and possible - we leave as an inquiry for the future. Occasionally 
note that though Mendez's Stirling numbers of the first and second kind 
are not Whitney numbers of an appropriate poset they do bear a striking 
resemblance to the latter. 

Whitney numbers, statistic, interpretation. It is well known [46] that 
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denoting set of n elements partition lattice by II„, the arrays 
are identified (see also Theorem 1.3 in [18]) as follows 



k 
r 



and 



n 

n — k 



w k (U n ) and j ™ } = W k (U n ) 



where tUfc(II n ) and PFjt(II n ) denote Whitney numbers of the first and sec- 
ond kind correspondingly. In order to recognize the possible evolvement of 
state of affairs while the combinatorics is concerned let us come back for a 
while to ^-extensions realm. There are several available ways to define com- 
binatorially |/«| an d arrays. Most of these ways are based on on 

static on set partitions (see for example [47], [19], [18], [14], [16], [17], [48], 
[49]). For example Gessel in [48] gave to |&| combinatorial interpretation 
as generating functions for an inversion statistics. In another source paper 
[49] Milne demonstrated that | | may be viewed in terms of inversions 
on partitions and that they count restricted growth functions using various 
statistics (see also [16]). We owe to Milne also the interpretation of |^| 
as sequences of lines in a corresponding vector space over finite field. In [19] 
Sagan delivered the major index statistics' interpretation of | ^ | array of 

^-Stirling numbers of the second kind. After that the authors of [18] con- 
structed a family p n (q) of posets as g-analogues of the set partition lattice 
(different from Dowling (/-analogue) in such a manner that (Theorem 5.3 in 
[18]) 



n 

n — k 



w k (p) and { } =W k (p) Vp £ p n {q) 

In — kin 



become Whithey numbers w k (p) and W k (p) of the first and second kind 
respectively. Whitney numbers for any graded poset may be looked at as 
Stirling like numbers. We shall indicate at the end of this survey a class of 
substantially new examples of such Stirling like arrays - after we inform on 
prefabs' structures. Meanwhile let us come back to the main challenge of 
^-extensions where we are faced with an ispiring obstacle. 
Surprise ? In pQ, [2] a family of the so called ^-Poisson processes was 
introduced i.e. the corresponding choice of the function sequence ip leads 
to the Poisson-like ^-Poisson process. Accordingly one would expect the 
extension of Dobinski formula to the ip- case - to be automatic. Of course 
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it makes no problem to call the numbers 

B n (tp) = ip n (x = l,ip) 

the ip - Bell numbers - whenever it makes sense - for example either the se- 
quence of these numbers has combinatorial interpretation and/or the defin- 
ing series below are convergent: 

fan 

B n {ip) = ip n (x = l,ip) = exp^ 1 TTf- (35) 

o<fc ^ 

The above might be a far reaching generalization of the standard case [23]. 
For example - what about the spectacularly natural and number theoretic 
important choice: ip = (y^) n >o {{F n ) - Fibonacci sequence )? In this con- 
nection (Fibonacci binomial coefficients [50] are natural numbers!) a ques- 
tion arises whether one can prescribe eventual arithmetic properties of some 
of "0-Beil numbers beyond g-extensions to any kind of composite modules 
as in [24], [51] or [52] and [53, 54] - see references therein. The papers just 
mentioned perform their investigation mostly in umbra symbolic Blissard 
language (see the introduction in [5] for historical remarks on Blissard's 
calculus roots). Note then (see On ^-extension remark above) that in 
the ^-extensions realm one may formally introduce the '(/'-extended umbra 
symbol B^ by analogy to the Bell's source of the idea article [41] as follows 

£(Y>)n+i = (£(Y0 +^ If, ^ = x k . (tP-B -umbra) 

(see symbolic formula (4.7) in [41] p. 264). The above definition is equivalent 
to 

BWWi = £ (f\ BW*> n>0. 

For the meaning of the ^-shift "+^" operator - already implicit in Ward's 
paper [55] - see [1, 2, 26, 31, 33]. See occasionally substantial reference to 
Ward [55] in Wagner's article [15] on generalized Stirling and Lah numbers. 
Question. Summarizing the discussion above - would we then - beyond the 
g-umbral extensions' realm - would we have ^-Bell numbers with Poisson - 
like processes background - and not related to a kind of Stirling numbers ex- 
tension - at least in a way we are acquainted with? Or should we introduce 
extended Stirling numbers in another way so as to be not related to Poisson 
- like processes beyond the g-umbral extensions' realm? 
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On this occasion note also that all ip extensions of umbral calculus do not ex- 
haust all possible representations of GHW. For GHW (Graves-Heisenberg- 
Weyl) algebra the most general representation of its' defining commutation 
relation is already implicit in [56] which serves [1, 2, 33] as the algebraic 
operator formulation of Ward's calculus of sequences [55]. Namely from the 
Rodrigues formula (Theorem 4.3. in [56]) with 

xq n „i{x) = q n (x) , Qq n = nq n -i 

it follows that 

[Q,x] = 1 , x = xQ'^ 1 

where Q - a differential operator [56] is a linear operator lowering degree of 

any polynomial by one. Q needs not to be a delta neither ^-delta operator 
[2, 1, 33]. We deal with such a case after the choice of admissible sequence 
[56, 1, 2, 33] different from ^-sequence (^y)ri>o or (^jr)n>o m an d 
(L~) for (_B~(^/>) Bell numbers. Then from [57] we know that basis con- 
sisting of the persistent root polynomials (ipk(x))k>o does not correspond 
to V-delta operator. However it determines [56] a differential operator i.e. 
the linear one lowering degree of any polynomial by one. Another possible 
"rescue" in seeking for the convenient, efficient structure with natural ob- 
jects corresponding to Stirling or Bell numbers and Dobinski-like formulas in 
special cases are the exponential structures and prefabs. For example read- 
ing [30] one notices (section 3) that the g-analog of the Stirling numbers of 
the second kind description developed by Morrison (compare with Section 
4 in [15] to see in which way it is complementary) constitutes the same 
example of Ward'ian - prefab'ian" extension as the Bender - Goldman [58] 
prefab example to be considered next right now. As noticed by Morrison 
the relevant prefab exponential formula may equally well be derived from 
the corresponding Stanley's exponential formula in [59] . 
Exponential structures versus prefabs. A subcase of Two General 
Classes. Exponential structures and exponential prefabs are - in Stanley's 
words - basically two ways of looking at the same phenomenon [59] . Before 
coming over to inspect [59] from the "Stirling point of view" let us give at 
first a family of decisive examples showing that prefabs are all around us in 
combinatorics especially when quite free extensions of Stirling numbers are 
concerned. The following example contains such a family. 
Bender - Goldman - Wagner Ward - prefab example. If corre- 
sponding "prefabian" g-Bell numbers B^[ e ^ (7) are defined as sums over k of 
S q (n,k) Stirling numbers of the (/-lattice of unordered direct sums decom- 
positions of the n-dimensional vector space V q ^ n over GF(q) = F q in sect. 
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2 of [15] then the formula (2.5) in [15] shows up equivalent to the Bender- 
Goldman exponential formula (17) from [58] - the source paper on prefabs 
- and in our ■(/'-extensions' notation formula (17) with D n (q) from [58] now 
reads: 

B?*f(x) = £ B^ffr) — = exp{exp,(x) - 1}. ( 7 - e.g.f.) 

n>0 1 

Here 

n 7 ! = (q n - l)(q n - q l )...{q n - q n ~ l ) = \GL n (F q )\, 

D Q (q) = 1 by convention while D n (q) = Bn Te ^ ( 7 ) = number of unordered 
direct sums decompositions of the vector space V q . n . Compare with formulae 
(4-6) in [60] representing the completely new class of combinatorial prefab 
structures with noncommutative and nonassociative composition (synthe- 
sis) of its objects. Note "The natural hint" on ip extensions remark there 
right below these formulae. Coming back to the Bender - Goldman - Wag- 
ner Ward - prefab example it is to be noticed that this is a special case of 
the First Class formula according to the terminology of three paths of gen- 
eralizations being developed in [15]. According to us Wagner justly refers 
his First class to Ward [55]. We propose to call this Wagner's First Class a 
"Ward'ian - prefab'ian" Class of extensions as the characterization formula 
(1.15) in [15] after being summed over k yields exactly ip- extension [60] 
of prefab exponential formula (12) from [58] where = (j-) n >o m Wag- 
ner's notation [15]. Note that our notation [1, 2, 33] is consequently always 
"Ward'ian". We also advocate by means of the present paper the attitude of 
Two General Classes. The Wagner's Class I is in our terms "Ward'ian - 
prefab'ian" (see (7 — e.g.f.) above) with F(n, k) Ward-prefab Stirling num- 
bers and with S q (n,k) as example. The second general class in our terms 
is "Newtonian" and it incorporates Wagner's Class II and Class III with 
N-W-C Stirling numbers U(n,k) and with Newtonian S q {n,k) and Gould- 
Carlitz-Milne S q (n,k) as examples mutually expressible each by the other 
one with help of re-scaling. One may see that really we are dealing here 
with the Newtonian way notifying that our N-W-C formula (20) extends 
(1.12) from the Class III of [15] and our N-W-C formula (14) extends (1.10) 
from the Class III of [15]. Note also that (5.28) from [15] via re-scaling 
coincides with q - N-W-C case of (20) i.e. with (q — Dobinski) formula. 
Here inevitable questions arise. For the Newtonian General class we have 
the extension (20) of Dobinski formula. So what about the corespondent 
formula for the Ward'ian - prefab'ian general class? ... And what about 
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The Two General Ways of this paper? The one "easy" way is Newto- 
nian. The other way seams to contain the g-extension as a kind of "singular 
point" in the domain of ?/;-umbral extensions. Is there at and beyond this 
"singular point" of the second way an another non-Newtonian second path 

- all-embracing what was left beyond the first way ? Before an attempt to 
answer some of these questions let us encourage ourselves by just recalling 
another distinguished example. 

This another crucial "Ward'ian - prefab'ian" example we owe to Gessel [48] 
with his (/-analog of the exponential formula as expressed by the Theorem 
5.2 from [48]. 

We also recall that the g-analog of the Stirling numbers of the second kind 
investigated by Morrison in Section 3 of [30] constitute the same example of 
Ward'ian - prefab'ian extension as in the Bender - Goldman - Wagner Ward 

- prefab example. As noticed there by Morrison the (7 — e.g.f.) prefab 
exponential formula may equally well be derived from the corresponding 
Stanley's exponential formula in [59]. Let us then now come over to these 
exponential structures of Stanley with an expected impact on the current 
considerations ( for definitions, theorems etc. see [59]). In this connection we 
recall quoting (notation from [59]) an important class of Stanley's Stirling - 
like numbers M J^ of the second and those of the first kind Stanley's Stirling 

- like numbers ■ Both kinds are characteristic immanent for counting of 
exponential structures (or equivalently - corresponding exponential prefabs) 
and inheriting from there their combinatorial meaning. This is due to the 
fact [59] that "with each exponential structure is associated an "exponential 
formula" and more generally a "convolution formula" which is an analogue 
of the well known exponential formula of enumerative combinatorics" [59]. 
Consequently with each exponential structure are associated Stirling-like , 
Bell-like numbers and Dobinski - like formulas are expected also. 
Exponential structures. Let {Q n }n>o be any exponential structure and 
let {M(n)} n >o be its denominator sequence i.e. M(n) = number of minimal 
elements of Q n . Let |Q n | be the number of elements of the poset Q n 

\Qn\ = E L 

Example: For Q = (n n ) n >i where II ra is the partition lattice of [n] we have 
M(n) = 1. 

Define "Whitney-Stanley" number S n k to be the number of ir £ Q n of 
degree equal to k > 1 i.e. 
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s n ,k — i- 

7r6Q n ,|7r|=fc 

Define S n> k - generating characteristic polynomials (vide exponential poly- 
nomials) in standard way 

n 

W n ( X )= X^=Y,Sn,kX k . 
irGQn k=l 

Then the exponential formula (Wq(x) = 1 = M(0)) becomes 
^ W n (x)y n _! 

n=0 v ; 

where 

n=l 

with the obvious identification of ^-extension choice here. Hence the poly- 
nomial sequence {p n (x) = )w>o constitutes the sequence of binomial 
polynomials i.e. the basic sequence of the corresponding delta operator 
Q = q(D). We observe then that 

n q k n 
Pn(x)=£^=£[0,l,2,...,Mn]* fe 
k=0 ^ > k=0 

are just exponential polynomials' sequence for the equidistant nodes case i.e. 
Newton-Stirling numbers of the second kind 5~ fc = ■ Both numbers 
and the exponential sequence are being bi-univocally determined by the 
exponential structure Q. This is a special case of the one already considered 
and we have as in this "Lupas case" the Newton-Stirling-Dobinski formula: 

Pn(x) —J] = £[0, 1, 2, k; b n ]x\ (N-S-Dob) 

exp W k =o K - k=0 

where (b n ) n >o is defined by 

n 
k=0 
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Note the identification b n (x) = jjj^, where 

Wn(x) = ~ M0,Vr)AW. 

\x is Mobius function and is unique minimal element adjoined to Q n . 
Corresponding Bell-like numbers are then given by 

1 00 u (t,\ n 

Pn{1) = Eir = B 011 ' 2 ' k] M ' {N ~ S - Bell) - 

exp(xj fc=o fc. fc=Q 

Besides those above - in Stanley's paper [59] there are implicitly present also 
inverse-dual "Whitney-Stanley" numbers s n ^ of the first kind i.e. 

Sn,k = ~ ^ M(0i 7r )- 

TreQ n ,\Tr\=k 

On this occasion and to the end of considerations on exponential structures 
and Stirling like numbers let us make few remarks, (/-extension of exponen- 
tial formula applied to enumeration of permutations by inversions is to be 
find in Gessel's paper [48] (see there Theorem 5.2.) where among others he 
naturally arrives at the (/-Stirling numbers of the first kind giving to them 
combinatorial interpretation. Recent extensions of the exponential formula 
in the prefab language [58] are to be find in [60]. Then note: exponen- 
tial structures, prefab exponential structures (extended ones - included) i.e. 
schemas where exponential formula holds-imply the existence of Stirling like 
and Bell like numbers. As for the Dobinski-like formulas one needs bino- 
mial or extended binomial coefficients' convolution as it is the case with 
-^-extensions of umbral calculus in its operator form. 

Information. On the basis of [60] the present author introduces new prefab 
posets' Whitney numbers in [61]. Two extreme in a sense constructions are 
proposed there. Namely the author of [61] introduced two natural partial 
orders: one < in grading-natural subsets of cobweb's prefabs sets [60] and 
in the second proposal one endows the set sums of the so called "prefabi- 
ants" with such another partial order that one arrives at Bell-like numbers 
including Fibonacci triad sequences introduced by the present author in [62]. 

A.3 

Other Generalizations in brief. We indicate here three kinds of exten- 
sions of Stirling and Bell numbers - including those which appear in coherent 
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states' applications in quantum optics on one side or in the extended rook 
theory on the other side. In the supplement for this brief account to follow 
on this topics let us note that apart from applications to extended coherent 
states' physics of quantum oscillators or strings [6 - 11, 63, 64] and related 
Feymann diagrams' description [65] where we face the spectacular and in- 
evitable emergence of extended Stirling and Bell numbers (consult also [66]) 
there exists a good deal of work done on discretization of space - time [67] 
and/or Schrodinger equation using umbral methods [68] and GHW algebra 
representations in particular (see: [67, 68] for references). 

A. 3. 1. An analog of logarithmic algebra. An extension of binomiality 
property from the algebra of formal power series to the algebra of formal 
Laurent series and then beyond leading to the Loeb's [69] iterated loga- 
rithmic algebra - was realized by Roman [70, 71] with the basic Logarithmic 
Binomial Formula at the start. The Logarithmic Binomial Formula and the 
iterated logarithmic algebra may be given their ^-analog including q-analog 
of the Logarithmic Fib-binomial Formula as shown in [72]. 

The extension of the iterated logarithmic algebra from [69] is the loga- 
rithmic algebra of Loeb-Rota [73] . This generalization of the formal Laurent 
series algebra retains main features and structure of an umbral calculus. 
Among others it allows for logarithmic analog of Appell or Sheffer polyno- 
mials. These give rise to Stirling- type formulas already in [69]. In [74] 
Kholodov has invented an analog of the logarithmic algebra from [69] in the 
shape of an umbral calculus on logarithmic algebras. Specifically (Example 
3.1 in [74]) the basic logarithmic algebra constructed via Jackson derivative 
d q gives rise to the analog of (/-Stirling formula. The mnemonic natural 
question arises: are similar constructs performable for ^-derivatives ? 

A. 3. II. Milne's Dobinski formula. In the classical umbral calculus rep- 
resented by the finite operator calculus of Rota the clue and source example 
of delta operator is 

A = exp{D} - I. 

Naturally such delta operator generates Stirling numbers of the second kind 
via 

Accordingly clue and source example of delta operator of the ip~ umbral 
calculus would be 

= exp{d^,} - I. 
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However already Milne's (/-extension [12] - contrary to the above - does not 
rely on Jackson derivative d q and it reads 




\k n\ 



where A k - th k-th difference operator is defined inductively so that 



A k q = (exp{D}-q k - 1 I)(exp{D}-q k - 2 I)(exp{D}-q k - 3 I)...( y exp{D}-q I) 



The corresponding (/-Dobinski formula ( 1.26 in [12]) looks " ^-familiar' 
(see: (35) above): 



The obvious challenge is an eventual application of that type extension to 
other umbral calculi - including the analog of the logarithmic algebra from 
[73] with use of < \ \ perhaps. 



A. 3. III. Normal ordering accomplishment and generating func- 
tions as coherent states. While staying with formal power series algebra 
or even its subalgebra of polynomials - still valuable extension have been 
applied as desired from at least two points of view: a) statistics , b) normal 
ordering task for quantum oscillator and strings. 

a) As for statistics recall that (/-Stirling numbers of the second kind may be 
treated as generating functions for various statistics counting ( see: [15] and 
references therein) . This type of role has been given by Wachs and White in 
[16] to p, q- Stirling number S Ptq (n,k) which is generating function for the 
two different joint distribution set partitions statistics. Wachs and White 
in [16] also hade proposed interpretations of their p, (/-analogue of Stirling 
numbers in terms of rook placements and restricted growth functions. From 
the defining recurrence (4) in [16] one sees that S\ >q (n,k) = < n k \ . 



b) Similar two parameter r, s-Stirling numbers S r ^ s (n, k\q) arise in the nor- 
mal ordering accomplishment for the expression [(a + ) r a s ] n [6,7] where a + 
and a stay for creation and annihilation operators for (/-deformed quan- 
tum oscillator which is equivalent to say that aa + — qa + a = 1. (For ex- 
ample a = d q , a + = x). From the recurrence (50) in [7] one sees that 

S rt i(n, k\q) = I n k \ . This special case of Gould-Carlitz-Milne (/-Stirling 



numbers < n k \ appears in [8]. The method to use it in order to recognize 




(36) 
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coherent states as combinatorial objects was invented by Katriel in [6]. The 
authors of [9-11, 63, 64] develop the consistent scheme of applications of the 
properties of S rtS (n, k\q = 1) = S r ^ s (n, k). These include [11] closed-form ex- 
pressions for S r;S (n,k), recursion relations, generating functions, Dobinski- 
type formulas. Recall that generating functions are identified with special 
expectation values in boson coherent states. Recall that S r ^ s (n,k) in ter- 
minology proposed in this paper are the special case of Newton-Stirling 
numbers of the second kind and correspondingly - Dobinski-type formulas 
are - see the Coherent States' Example in A. 2. above. A new perspective 
opens while considering normal ordering task not only for quantum oscilla- 
tor but also for strings which are " many, many oscillators" . The first steps 
had been spectacular accomplished by the authors of [63]. These authors 
obtained not only analytical expressions but also a combinatorial interpreta- 
tion of the corresponding "very much extended" Stirling and Bell numbers. 
Their properties are interpreted in [63] in terms of specific graphs. At the 
same time the authors of [63] consider an invention of a g-analog of [63] to 
be "an outstanding problem". 

A. 3. IV. From Howard via Hsu and Yu and Shiue to Remmenl and 
Wachs extensions. Information in brief. Howard's [75], via Hsu's and 
Yu's [76] and Shiue's [77] to Remmenl and Wachs [78] sequence of ex- 
tensions starts with degenerate weighted Stirling numbers [75] used later 
on by the authors of [76, 77] to propose respectable, unified approach to 
generalized Stirling numbers. This sequence of extensions ends with elab- 
orated extended rook theory [78] with its generalized Stirling numbers and 
(p, g)-analogues of Hsu and Shiue extensions. Recall that (p, g)-analogues 
of Stirling numbers were introduced by Wachs and White in [16]. In more 
detail. Hsu and Shiue had provided a unified scheme for many extensions 
of Stirling numbers of both kinds known before [77]. They introduced cor- 
responding unified extensions under the notation: 

\ 2 

Sn,k{ a iPi r ) and S n,ki a 'P' r )' 

and ^,(1,0,0) = {"} andS^KAr) = 
5'n,fc(/3>a>- r )- 

Guided by Wachs and White ideas from [16] Remmel and Wachs have defined 
in [78] two natural p, (/-analogues of Hsu and Shiue extensions 

S^,fc(a,/V), * = 1,2. 



such that S U:k ( 1,0,0) 



n 
k 
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For that to do Remmel and Wachs have used what we would call the ip Ptq 
admissible sequence ip Pjq = ([n] Pi? ) n >o where [y] Ptq = P p Z q which for 7 € 
N U {0} becomes the known [16] extension of Gauss extension i.e. [n] Ptq = 
qU-i _|_ _|_ ^ _|_pn-2g _|_^n-i_ F ac torials and ^g-binomial coefficients 

are then defined accordingly naturally (see: A.l. Notation.) 

A. 4. Extended umbral calculus and some corresponding exten- 
sions of Stirling and Bell numbers. Further examples. In this part 
of our presentation we just list some examples at hand where evidently ip- 
extension is behind the scenario for special ^-admissible sequence choices. 
Example A. 4.1. In Katriel and Kibler's celebrated paper [8] on normal or- 
dering for deformed boson operators and operator- valued deformed Stirling 
numbers one uses the following ^-admissible sequence 

1 Q n _ v n 

( «■ = w >«.. H = Vr^. 

from Wachs and White source paper [16]. 

Example A. 4. 2. In Schork's paper [79] on fermionic relatives of Stirling 
and Lah numbers one uses the following ^-admissible sequence 

1 1 — (— l) n 
( i>n = |^pj)n>0, [nf = n q= -i = . 

Example A. 4. 3. In Parthasarathy's paper [80] on fermionic numbers and 
their roles in some physical problems one uses the following V-admissible 
sequence 



1 x _ l-(-l) n g n 



( 4>n = TT7t)™>0> i n 



[n]/! /n " U ' LJ l + q ' 

Here in [80] the g-fermion numbers emerging from the g-fermion oscillator 
algebra are used to reproduce the g-fermionic Stirling and Bell numbers. 
New recurrence relations for the expansion coefficients in the 'anti-normal 
ordering' of the q-fermion operators are derived. Corresponding ^-extended 
Dobinski formula (see: (15) in [80]) is derived. 

Example A. 4. 4. In the paper [11] on extended Bell and Stirling numbers 
from hypergeometric exponentiation one uses the following ^-admissible se- 
quence 

Here in [11] elements of ip L -umbral calculus are at work. Among others the 
corresponding ^-extended Dobinski formula (see: (15) in [11]) is derived. 
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Example A. 4. 5. In the paper [81] on representations of the so called 
" Monomiality Principle" with Sheffer-type polynomials and boson normal 
ordering - just the standard ^-admissible sequence choice 

( V>n = A)n>0- 

of the classical non-extended umbral calculus is naturally abounding in un- 
countable formal series indicators of delta operators examples; see: a) - 
g) page 3 in [81]. Note there also GHW - algebra formula (12). As for 
what the authors re-discover (?) to be the so called "monomiality principle" 
one should note and compare this with the source paper [56] from 1978 by 
George Markowsky on " Differential operators and Theory of Binomial Enu- 
meration". In particular see the GHW - algebra in spirit Theorem 4.3 in 
[56] (see also [1, 33] for more on that). 

A. 5. Historical and bibliographical relevant remarks. To this end we 
shall here list few peculiar relevant remarks of historical and bibliographical 
character. 

A. 5.1. Remark The history of GHW algebra has its roots not later then 
since Graves' work [32] "On the principles which regulate the interchange 
of symbols in certain symbolic equations" from (1853-1857). See [82] by 
O.V. Viskov "On One Result of George Boole" (in Russian) from 1997. 
A. 5. 2. Remark Generalizations given by formulas (3) and (4) from Cakic 
and Milovanovic paper [83] for another extensions of Stirling numbers of the 
second kind as well as their related properties are an old result published 
in d'Ocagne M. article in 1887 [84]. Many other later generalizations (see 
[83]) are consequence of Chak's work [85] and Toscano papers [86-88]. The 
relevant papers of importance (see: [83]) are those [89-93] and [75]. 
A. 5. 3. Remark The relevance of Schlomilch's work [94] from 1852 to N-W- 
C Stirling numbers is taken down here with pleasure. Another interesting 
paper refereeing directly to the original Dobinski's work [4] and Dobinski's 
point of view is the Fekete's paper [95] from 1999. 
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